Energy Systematics of Jet Tomography at RHIC: y/s = 62.4 vs 200 AGeV 
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The collision energy dependence of jet tomography is investigated within the GLV formalism. 
The emphasis is on estimating systematic uncertainties resulting from the interplay of energy loss 
fluctuations and the rapid increase of the parton transverse momentum slopes as y/s decreases from 
200 to 62 AGeV. 
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I. INTRODUCTION 

We study the energy dependence systematics of 
jet tomography in nuclear collisions within the GLV 
formalism|l Q in the y/s = 62 to 200 AGeV range. In 
general, tomography provides a map of the density of 
matter through the study of the angular pattern of at- 
tenuation of a calibrated source beam. Extensive high 
Pt data on single and dijet azimuthal correlations and 
relative to the reaction plane have already been com- 
piled at 200 AGeV 0, and used to determine opacity 
and density of the quark-gluon plasma (QGP) produced 
in central Au+Au reactions^. The interest in applying 
the same techniques to lower 62 AGeV data |6J is to test 
experimentally the predicted y/s decrease of the QGP 
density and the predicted variation of the gluon/quark 
jet source. In this paper, we calculate the nuclear modi- 
fication factor, Raa(pt, V = 0, y/s), for central Au — Au 
collisions at both y/s = 62 A, 200 AGeV for neutral pions. 
Previous predictions for 62 AGeV have been published by 
Wang0 and Vitev@. 

We concentrate here on the role of energy loss fluc- 
tuations 0, El on the predicted single hadron attenu- 
ation pattern in order to gain an estimate of some of 
the systematic theoretical errors in the jet tomographic 
technique. To isolate the role of fluctuations we neglect 
kx smearing, the Cronin enhancement, gluon and quark 
shadowing, and nonperturbative baryon dynamical con- 
tributions that strongly distort the hadron spectra be- 
low px < 4 — 5 GeVjj]. Therefore, our calculations can 
only be compared to unidentified charged particles out- 
side the baryon anomaly region^, i.e. pt ^ 5 GeV. Even 
for the identified ir° spectra, the calculations neglecting 
the rapidly rising Cronin effect with decreasing energy 
limits the applicability to the pr > 4 GeV region. See 
different estimates of those effects in Refs.0,U The in- 
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termediate pr is also strongly influenced by the trans- 
verse Dopplcr shifted hydrodynamic flow of the QGP 
bulk matter |lfj as well as possible coalescence hadroniza- 
tion mechanisms |l lj not described by attenuated parton 
fragmentation. 



The tomographic information about the maximum 
densities of the QGP attained at RHIC as a function of 
energy is best isolated from the pr > 5 GeV region. How- 
ever, it is important to examine critically the underlying 
assumptions of jet tomography to assess the theoretical 
systematic errors in its application to nuclear reactions. 
We test below the influence of the shape of the energy 
loss fraction spectrum, P(e, e) about the mean energy 
loss fraction, e. The tomographic information is encoded 
in e because for Bjorken expansion it is is proportional to 
the produced gluon rapidity density, dN g /dy as well as 
the path length. Neglecting fluctuations as in |l2|. i.e., 
P(e, e) = 5(e — e), a fit of Raa provides a biased measure 
of e and hence dN g /dy. 



As first emphasized in [13J fluctuations about e must 
be considered in realistic applications because the very 
rapid decrease of the transverse momentum spectra of 
partons induces a bias to lower effective e. In Ref.0 we 
quantified the magnitude of the fluctuation bias in the 
GLV approach through a renormalization factor Z . We 
showed that for 130 AGeV conditions Z ~ 0.5 due to this 
bias, and hence the fitted e from Raa must be upward 
corrected by a factor 1/Z s» 2. 



At lower energies, the high px slopes of the transverse 
momentum spectra increase more and more as y/s de- 
creases as seen in Fig. 1. Thus the calibrated gluon and 
quark source beams change systematically as a function 
of y/s. This increases the fluctuation biases and hence 
it is important to investigate how sensitive is the tomo- 
graphic analysis to not only e but also to the shape of 
that distribution. 
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FIG. 1: Invariant cross-sections for tv° production in p-p 
collisions at ^/s = 30.7,53.1,62,4,200 GeV. Data is from 
PHENIX0 (200 GeV) and ISR (all other energies) collab- 
orations. 



III. FINAL STATE QUENCHING AND 
FLUCTUATIONS 

In a dense QCD medium the induced radiative energy 
loss reduces the initial transverse momentum p T of a pro- 
duced hard parton prior to hadron fragmentation. The 
reduction factor (1 — e) depends in general on both the 
initial parton rapidity y and transverse momentum p T 
as well as on local QGP density that varies with s and 
impact parameter b 

In this paper, we consider only mid-rapidity (y = 0) 
and central b << R collisions and further simplify the 
calculation by neglecting the slow px dependence of the 
energy loss fraction, e. The parton before fragmentation 
has pip — p T (l — e) in this picture. We neglect acopla- 
narity as well. 

Fluctuation in the energy loss is taken into account via 
a probability function P(e, e) using 0] 

Eh ^T = K ^ j dXldX2 f-/A(^Q 2 )h/A(x2,Q 2 ) 

dt J z c nz c 



II. CALCULATION OF SPECTRA 

The neutral pion invariant inclusive cross section in pp 
at various values of a/s can be calculated via conventional 
collinear factorized pQCD, 

f 

Eh ~dFp Jdxidx 2 f a /A(xi,Q 2 )fb/A{x2,Q 2 ) 

da ab^cd j^ 0/e ( Ze; Q2) 

dt irz c 

The inclusive number distribution, dN 7 * /dyd 2 pT, in 
A + A collisions is obtained in the absence of nuclear ef- 
fects by simply multiplying Eq. by the Glauber ge- 
ometric (binary collision density) factor, TAB{b)- We 
use standard Leading Order MRS D- parton distribu- 
tion functions f a (x) (as in Jl7ll ). and KKP fragmentation 
functions Dh(z c ) (as in 18]). Here, x± is the momentum 
fraction of the projectile parton in the collision while x 2 is 
the momentum fraction of the target parton. The K (-Js) 
is a factor that simulates higher order effects. 

Figure compares the invariant pQCD cross-sections 
for 7T° production to available data. The particular values 
used for K and Q 2 for the fits are given K = 2.0, Q 2 = 
p 2 T for 200 GeV, K = 2.5, Q 2 = 0.5p| for 62.4 GeV, 
K = 3.5, Q 2 = 0.5p^ for 53.1 GeV and K = 4.5, Q 2 = 
O.^Pt for 30.7 GeV. The rapid increase of the slopes with 
decreasing t/s is well reproduced by pQCD. However, the 
normalization corrections due to NLO processes increases 
as s decreases. 



where z* — Pt-k/Ptc ~ z c/(l — e )- The parameter e de- 
termines the magnitude of the quenching and is expected 
decrease as the local density oc dN g /dy decreases at lower 
energies. The e is interpreted as the average px shift. 

In and @, -P(e) was computed assuming Poisson 
fluctuations of the radiated gluon number. This distri- 
bution is thus characterized by the mean number of radi- 
ated gluons, (N g ), as well as the mean energy loss e. An 
inherent problem of treating fluctuations via a Poisson 
process is that it leads to violations of the energy bound 
e = Y^n=i e ™ < 1- F° r small 1 g < 0.5 this is a small effect 
that can be approximately corrected by truncating the 
distribution at e = 1. However, for high opacity QGP 
matter produced at RHIC, the gluon jet e g ~ 0.6 — 0.8 
approaches unityp| and the truncation scheme leads to 
an approximately uniform distribution limiting e g < 0.5. 

The truncated Poisson ansatz for fluctuations provides 
a natural way to impose unitary and kinematic con- 
straints on the induced radiative energy loss. However, 
it is possible that other distributions could be operative 
in the opaque limit. In this paper we explore two simpli- 
fied forms of fluctuation distributions to assess some of 
the the systematic uncertainties in the predicted quench- 
ing factors. One is a simplified "uniform" model that 
reproduces the essential features of the truncated Pois- 
son ansatz of 0, @- The second distribution is called 
"squeezed" because it accumulates strength near e ~ 1 
in the opaque limit. 

The "uniform" distribution is assumed to have the 
form 

[ ' ' \ 1 if 0.5 <e 
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FIG. 2: Raa{pt,V* = 200) with squeezed fluctuations FIG. 3: Predicted energy evolution of Raa with squeezed 
The 7T° data are from PHENIX|| fluctuations at = 62.4, 200 GeV. 



The "squeezed" distribution is assumed to have the form 



8(0<e<2l) 
2e 

fl(2e-Kc<l) 
2(1-*) 



if < e < 0.5 
if 0.5 < e 



Multiple gluon convolutions lead to phase space suppres- 
sion e™ of low total energy loss, which is mimicked in 
© by a simple 9 function. A more detailed analysis 
of branching versus truncated Poisson distributions will 
be presented elsewhere [Tflj. For our present purposes of 
gaining a first measure of systematic errors the above 
overly simplified forms are adequate. Note that while 
the uniform distribution has a limiting 1 < \ even for 
gluons, the "squeezed" one can easily accomodate the 
more realistic i g > 0.5 situation at RHIC. 

The shape of P(e, e) depends not only on e but also on 

< N g >, the average number of total gluons emitted. A 
smaller value of < N g > weights the lower n values more 
and leads to a more uniform P(e, e) while larger values of 

< N g > weight higher n values [H| ■ The average number 
of gluons radiated increased with ^fs. Thus, we expect 
the uniform distribution to be a better approximation 
for P(e, e) at 62.4 GeV while the "squeezed" distribution 
could be more realistic at 200 GeV. 

The nuclear modification factor can then be calculated 
for either distribution from Eq. [21 via 



Raa(pt, s) 



da x o(e(s)) 
da^o (0) 



(3) 



A. Raa with Squeezed Fluctuations 



Our first step is to fit e(200) to the 200 AGeV PHENIX 
data[| in the region pr ~ 8 GeV. Figure [3 shows the re- 
sulting Raa- The quenched pQCD prediction for lower 



Pt < 6 GeV increasingly deviate from the data because 
intrinsic kx smearing and Cronin calculations are ne- 
glected here. However, the approximate pt independence 
of higher pr is well accounted for as in Q ■ The band cor- 
responds to an assumed average mean energy loss frac- 
tion for gluons at 200 AGeV 0.65 < e g (200) < 0.76 using 
the squeezed distribution. Further below we will contrast 
the results to the energy evolution with the uniform dis- 
tribution. 

The beam energy dependence of e(s) is predicted to 
follow the energy dependence of the produced dN g /dy. 
We also take into account the Casimir difference between 
quark and gluon energy loss. The mean energy loss frac- 
tion of parton type c = q,g evolves with energy as 



C g \dNg(200)/dy J 



(4) 



Thus e q = 4/9e g 

The predicted Raa at 62 AGeV is obtained using the 
multiplicity systematics from PHOBOS[lJ| that suggests 
a value for dN tot dy « 650-770. From entropy conserva- 
tion we expect that the gluon rapidity density is compa- 
rably decrease from about 1000 at 200 AGeV. Therefore, 
we expect that e 9 (62.4) » 0.65-0.77 e g (200). The bands 
of Raa predictions at 62.4 GeV are shown in Figure 
The ranges of average quenching shown in the figure are 
0.42 < e(62.4) < 0.49 and 0.50 < e(62.4) < 0.59 for mul- 
tiplicity values of 650 and 770 respectively. These ranges 
are calculated by realizing that ^^-(200) w 1000. 



B. Raa with Uniform Fluctuations 

We now repeat the Raa calculations in the previous 
section using the uniform distribution for P(e, e) from 
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FIG. 4: Raa at 200 GeV calculated using the uniform distri- FIG ' 5: Raa at V 3 = 62 ' 4 ' 200 GeV calc ulated using 
bution. The tt° data are from PHENIXH. unlform distribution. 



the 



Eq.@. This is done to investigate the sensitivity of the 
results to the form of this distribution. The uniform 
distribution is the same as the squeezed distribution for 
€ c < 0.5. Therefore the quarks energy loss fraction distri- 
bution is identical for both distributions. However, since 
the gluon GLV energy loss fraction e g > 0.5 in this en- 
ergy range, the uniform distribution actually saturates 
the mean energy loss fraction of gluons at e g — 1/2. This 
contrasts with the squeezed distribution that squeezes 
the gluon fraction energy loss toward the e = 1 bound- 
ary. Therefore the gluon quenching is reduced with the 
uniform distribution while preserving the quark quench- 
ing. As we note further below, this also changes the glue 
to quark contribution ratios at high pr- 

Raa is compared to PHENIX data using the uniform 
fluctuation distribution in Figure 0] The curves are dis- 
tinctly flatter than the ones in the previous section (Fig- 
ure This is because the uniform quenching distribu- 
tion does not suppress so much the gluon jets than the 
squeezed distribution. The band in this case corresponds 
to 0.70 < e g (200) < 0.80. In the uniform case this implies 
that the mean quark energy loss fraction at y/s = 200 is 
0.31 < e g (200) < 0.36 

Figurc[S]shows the Raa prediction for 62 AGeV in this 
case. The range of mean gluon energy fractional energy 
loss is 0.46 < e ff (62) < 0.52 and 0.54 < e g (62) < 0.62 
for the different estimates of the dN g /dy(62) = 650, 770 
at 62 AGeV respectively. The nuclear modification fac- 
tor has a weaker energy dependence in this case and is 
generally somewhat smaller relative to the squeezed fluc- 
tuation case discussed in the previous section. This is 
due to the approximate saturation property of e g < 1/2, 
as in the full Poisson fluctuation case [8j • 

These calculations for the nuclear modification factor 
at y/s — 62.4 GeV are consistent with the calculations 



performed in [8j. In the pt region 4-6 GeV the predic- 
tion bands of both the squeezed and uniform distribu- 
tions are within the approximate range as predicted by 
Vitev in |8( . The main differences in the predictions can 
be explained by noting that the calculation in Q includes 
models for Cronin interactions as well as initial parton fcy 
smearing which are not included in the current calcula- 
tion. Both of these effects tend to deplete hadrons in the 
higher transverse momentum regions and increase their 
number in the intermediate transverse momentum region 
from 2-4 GeV. This will explain the lower quenching at 
lower px predicted by Vitev as compared to the current 
calculation. 

Further comparisons can be attempted with recent 
data from the PHOBOS collaboration as seen in 6] . This 
publication, however, provides data at lower levels of px 
and measures unidentified charged hadron cross-sections. 
This leads to a double fold problem with the comparison. 
Firstly, our calculation does not apply in the kinematic 
region probed by PHOBOS in 6j. Lower values of pr 
are not modeled well by our calculation due to previ- 
ously mentioned reasons of neglecting Cronin effect and 
kx smearing. Also, the calculation of charged hadrons 
would only be similar to neutral pions where there were 
not many baryons present. The kinematic region probed 
by PHOBOS is dominated by baryonic data which is not 
taken into account in the calculation. This effect would 
increase Raa as can be seen from the data shown in ||| • 



IV. COMPARISON OF Raa FOR DIFFERENT 
VALUES OF y/s 

Note that the predicted Raa(pt, \A = 62) have a neg- 
ative px slope compared to the generally flat Raa at 200 
AGeV. At very high pr, the 62.4 GeV curves even dip be- 



low the ones at 200 GeV. This higher slope can be more 
robustly seen by calculating the variable R s (s) as 



"i 1 r 



R s (s) 



Raa(s) 
Raa(200) 



(5) 



showed in Figure El The figure shows multiple curves for 
both types of quench distributions. One can see that the 
R s curves have a distinct downward slope due to the in- 
creasing power law index of the underlying pQCD parton 
distributions at lower energies. These indices vary much 
more with p T at 62 AGeV than at 200 AGeV. Therefore, 
the downward slope with increasing pt is due to the in- 
creased slope of the input spectra with increasing pt at 
lower y/s. 

It is remarkable that even with a smaller average en- 
ergy loss at 62 AGeV, the more rapid decrease of high pt 
parton production compensates to make it appear that 
there is greater quenching. The y/s systematic varia- 
tion of the calibrated initial jet source transverse spectra 
used leads to this anti-intuitive feature and must be care- 
fully taken into account in future tomographic inversion 
of single and di-jet quaenching. Note that R s much less 
sensitive to the propagation of the uncertainty in Raa at 
200 AGeV because absolute dN g /dy uncertainty tends to 
cancel out in this R s ratio. This variable is also seen to 
be able to differentiate clearly between these two types 
of fluctuation distributions. 
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FIG. 6: The ratio of Raa at 62.4 GeV to the one at 200 GeV 
for different multiplicity assumptions. Curves for squeezed 
distributions tend to be higher than curves for uniform dis- 
tributions. 
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FIG. 7: The gluon and quark contributions to the total 
Raa{Pt,\/s — 200) are shown using a uniform distribution. 
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FIG. 8: The gluon and quark contributions to the total 
Raa(Pt, yfs = 62) are shown. In contrast to y/s = 200 AGeV 
the 7T° at 62 AGeV are dominated by quark fragmentation. 



V. QUARK VS GLUON FRAGMENTATION 

An interesting theoretical measure of the underlying 
dynamics is shown in fig. (|7l8|) 

We see that while quenched quark and gluon frag- 
mentation contributes comparable amounts to the final 
quenched w° spectrum at 200 AGeV, The gluon contri- 
bution is almost completely quenched by 62 AGeV in our 
calculations in spite of the smaller e g at 62. The quark 
and gluon p? transverse evolve with s in different cal- 
culable ways in pQCD. This evolution however causes 
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FIG. 9: The ratio of it from gluons both quenched and un- 
quenched at y/s = 62.4,200 GeV. ^-(62.4) is assumed to 
770. Lower levels of quenching are used to make this esti- 
mate. e S queezed is 0.65 and 0.50 at 200 GeV and 62.4 GeV 
respectively. e un iform is 0.7 and 0.54 at 200 GeV and 62.4 
GeV respectively. 



the rapid disappearance of the gluonic contribution as s 
decreases in the case of the uniform distribution. 

In FigEJwe see that the form of the fluctuation distri- 
bution strongly influences the pion from gluon jet frac- 
tion. For the squeezed distribution, the gluon jet contri- 
bution is already negligible at 200 AGeV and of course 
remains so at 62 modulo a small increase due to the de- 
crease of e(62) relative to 200. 



These results suggest that identified particle ratios sen- 
sitive to quark vs gluon fragmentation^^ in the high 
Pt > 5 GeV domain would provide valuable to constrain 
on the form of the energy loss fluctuations and to map 
out more accurately the interplay of quark and gluon 
contributions. 



VI. CONCLUSIONS 

Predictions for Raa {pt , y/s = 62) for 7r° were made 
with a focus on the systematic uncertainties resulting 
from as yet un-tested assumptions about the energy loss 
fluctuation spectrum. The most unintuitive aspect of the 
prediction is a decreasing Raa with px at lower energies 
that is due to the more rapid decrease of the produced 
parton px distributions. The negative pt slope of Raa 
provides an important constraint on the fluctuation dis- 
tributions. The pion fractions from gluon jets were also 
shown to be sensitive to the underlying fluctuation distri- 
butions. Identified particle ratios at pt > 5 GeV should 
help to further test the distribution. Such further studies 
are essential in order to calibrate more accurately the jet 
tomographic measure of the initial QGP densities pro- 
duced at RHIC. 
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